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Abstract: We consider the problem of unique identification of dielectric coefficients 
for gratings and sound speeds for wave guides from scattering data. We prove 
that the "propagating modes" given for all frequencies uniquely determine these 
coefficients. The gratings may contain conductors as well as dielectrics and the 
boundaries of the conductors are also determined by the propagating modes. 

§0. Introduction 

Consider Maxwell's equations for time-harmonic electric and magnetic fields, 
exp(— iu)t)E(xi, X2, X3) and exp(— iu>t)H(xx, £2, Xa), in the absence of currents and 
charges 

V x E — iu^qH = 0, 

V x H + iue(x)E = 0. 

In this paper we study the inverse problem of determining the electric permittivity, 
e, and hence the dielectric coefficient, e/eo, from scattering data for these equations. 
The fundamental assumptions are that e is independent of X3, 27r-periodic in x\ and 
constant (= e ) for \x%\ > T. These conditions are designed to model a dielectric 
"grating" extending throughout the region \x2\ < T (c.f. [P], [BDC]). We also allow 
for conducting bodies embedded in the dielectric as long as they satisfy conditions 
analogous to our conditions on e: their boundaries should be invariant with respect 
to all translations in xa and translations by 2ir in xx, and they should be contained 
in |a;2| < T. 

To define data sets for this inverse problem it is customary to consider the 
scattering problem for fields with either the transverse electric (TE) or transverse 
magnetic (TM) polarizations, respectively E(x) — (0, 0, u(xi, X2)) and H(x) = 
(0, 0, v(xi,X2))- These polarizations reduce Maxwell's equations, respectively, to 



'" Au + uj 2 [1qu = (TE) 



e(x) 
and 

V • 4tV« + w 2 Mo« = 0. (TM) 
e(x) 

In the case of embedded conductors we consider the TE polarization in the exterior 
of the conductors with the Dirichlct condition, u = 0, on the boundary, since this 
corresponds to E = and H ■ h = on the surfaces of the conductors. For our 
purposes it is convenient to write both (TE) and (TM) as 



Lu — k 2 u = 0, 
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where L = — A for \x 2 \ > T and k 2 = uj 2 ^ € o- 

We will also present the analogous inverse problem for acoustic wave guides. 
This requires only small modifications of the arguments for gratings. The wave 
guides that we consider are simply slabs, {0 < x\ < B}, in which the sound speed 
c is a function of (xi,^)- We assume that c(xi,x 2 ) — Co(xi) for \x 2 \ > T, and 
impose Dirichlet condition on x\ = 0, and the Neumann condition, d Xl u = 0, on 
x\ = B. These boundary conditions correspond to an acoustically soft reflecting 
surface at x\ = and an acoustically hard reflecting surface at x\ = B, modelling 
underwater sound propagation with x\ as depth (c.f. [BGWX]). We will show that 
scattering data from propagating modes for the operator L = — c 2 (x)A with these 
boundary conditions determine c(x). 

In both these settings we will apply recent results on inverse coefficient problems 
for hyperbolic equations (Belishev [B], Kachalov-Kurylev-Lassas [KKL] and Eskin 
[E1],[E2]). In those papers the data for the inverse problem is the Dirichlet-to- 
Neumann map. Hence the objective here will be to show that the scattering data 
determine the Dirichlct-to-Neumann map on a line x 2 = T . 

Inverse scattering problems for dielectric gratings have been studied previously 
in [BDC], [BF], [K], [HK] and [EY]. These articles consider primarily the inverse 
problem of finding the boundaries of conductors embedded in a dielectric of constant 
permitivity from scattering data. To the best of our knowledge the present paper 
is the first to show that a variable dielectric is uniquely determined by scattering 
data. 

Inverse coefficient problems for wave guides were studied in [BGMX], [GMX], 
[M] , [X] and [DM] . These papers give methods for recovering the sound speed. We 
only consider the uniqueness problem and prove that the sound speed is uniquely 
determined by the propagating modes. Our approach was influenced by the work 
of S.Dcdiu and J. McLaughlin, [DM], which also uses propagating modes. 

§1. Statement of results 

Our results for gratings hold under mild conditions on the operator on L in (1). 
We will assume that L is a second order elliptic operator on D C K 2 which is 
symmetric in the inner product 



The weight a(x) is continuous and strictly positive on D. The coefficients of L + A 
are supported in [a^l < T, and L commutes with translation by 2ir in x\. Likewise 
a(x) — 1 is supported in \x 2 \ < T and a{x\ + 2ir,x 2 ) — a(x\,x 2 ). We will also 
assume that the region D is invariant under translation by 27r in x\, and boundary 
D is smooth. There are two cases that we wish to consider. 

Case 1: D is connected and contains {\x 2 \ > T}. In other words, while there may 
be some holes in D, they do not disconnect D, and they are contained in \x 2 \ < T. 

Case 2: D is connected and we have the inclusions 




{x 2 > T} C D C {x 2 > -R} 



for some R > 0. 
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The domain of L will be H 2 (D)nH^ (D). By H*(D) (resp. H k (D)) we mean 
functions satisfying 

u( Xl +2tt,x 2 ) = e 2nia u{x u x 2 ) (1) 

such that <f>(x 1 )u(x 1 ,x 2 ) € H k (D) (resp. i? fc (L>)) for all <f> e C C °°(R). Note that 
this domain for L corresponds to the Dirichlet boundary condition on dD. 

For wave guides we simply take L = —c(x) 2 A on D = {x e M 2 : < x\ < 
B} 7 with Dirichlet and Neumann boundary conditions on ii =0 and X\ = B, 
respectively. As indicated above, c(x) = cq{x\) when \x 2 \ > T. 

For both gratings and wave guides scattering data at fixed energy k 2 are obtained 
from the "propagating modes". In the case of gratings the scattered wave v+ is 
obtained by solving (L — k 2 )u = in D with u = exp(Z • x) + v+, where I ■ I = k 2 
and v + is the "outgoing" solution to 

(L-k 2 )v = -(L-k 2 )e a - x 

obtained as the limit as Im{/c} — > 0+ (see below). For gratings we will only use 
scattering data from incident waves exp(il ■ x) which satisfy the condition (2), i.e. 
those with l\ = n + a,n e Z, l 2 = —\]k 2 — (m + a) 2 . For x 2 > T the scattered 
wave v+ has the form 

v+(x, n,k) = J2 e^ m+ ^ Xl+X2 V k2 -( m + a )^a m (n, k) (2) 
= S m+a)xi+X2 ^ k2 - {m+a)2] a m {n, k) + 0(e- Sx2 ) 

(m+a) 2 <fc 2 

for some 6 > 0, provided that k does not belong to the set of "thresholds ", {k : 
k 2 = (p + a) 2 ,p e Z}. We call {a m (n,k), n,meZ : (m + a) 2 < k 2 , (n + a) 2 < k 2 } 
the scattering data at energy k 2 from "propagating modes" . 

For wave guides, since we are taking L = — c (xi) 2 A as the unperturbed 
operator, the scattered wave v + (x,k,n) is obtained by solving (L — k 2 )u = 
with u = $(x,k,n) + v + , where $ is a generalized eigenfunction for L 0} i.e. 
&(x,k,ri) = cxp(-ix 2 y/ Hn{k))<j)n{xi, k), where 

< + k 2 c^ 2 cj, n = fi n (k)<f> n , 0„(O) - 0, <f/ n {B) = 0, and fi n (k) > 0. 

For x 2 > T the scattered wave v + has the form 

v+(x,k,n) = 6 m (fc,n)e^V^W0 m (x!,fc) 

b m (k,n)e^^ k ~Um(x 1 ,k) + 0(e- Sx2 ) (2') 

/j m (fc)>0 

for some S > 0, provided that k does not belong to the set of thresholds, {k : 
Hp(k) = 0,p G N}. Here we call {b m (k,n), n, m e N : fi n (k) > 0, /ti TO (fc) > 0} the 
scattering data at energy k 2 from propagating modes. 

Note that in these definitions the functions \Jk 2 — (m + a) 2 and \J \J, m {k) will be 
chosen so that they extend into Im{fc} > with positive imaginary parts. Letting 
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z(k) stand for either k 2 — (m + a) 2 or /i m (/c), this choice amounts to choosing 
yj z(k) > when z(k) > and k > 0, \/ < when z(fc) > and k < 0, and 
^/ z(/c) = iy/\z(k)\ when z(fc) < 0. We will follow these conventions in the rest of 
the paper. 

With the preceding definitions we have: 

Theorem 1: The scattering data from propagating modes in x 2 > T given for all 
k determine D and e(x) for gratings with either the (TE) or (TM) polarizations, 
and c(x) for wave guides. 

The proof of theorem will proceed as follows. We will consider "generalized 
distorted plane waves" 



which arc defined without the restrictions k 2 > (n + a) 2 and fi n (k) > 0. These 
generalized distorted plane waves exist for k € R\S, where S is a discrete set. Note 
that when k 2 < (n + a) 2 or yU„(fc) < these generalized distorted plane waves grow 
exponentially as x 2 — > oo. In §2 and §3 we show that the set of generalized distorted 
plane waves, given for a fixed k and all n, uniquely determines the Dirichlet-to- 
Neumann map on a suitable line x 2 — T for all choices of k outside a discrete 
set. We also show that, making use of the analytic continuation to Im{fc} > 
of the v + (x, k, n)'s, these generalized distorted plane waves are determined by the 
scattering data from propagating modes. Thus, under the hypotheses of Theorem 1, 
the Dirichlet-to-Neumann map is known on x 2 = T for all k £ R outside a discrete 
set. Since this is equivalent to knowing the hyperbolic Dirichlct-to-Neumann map 
for the wave equations utt = Lu on x-i = T, the proof of Theorem 1 will be reduced 
to the results on hyperbolic inverse coefficient problems cited above. Since analytic 
continuation plays a big role here, there are many variations on the set of k for 
which the propagating modes are known which lead to the same results. 

§2. Determination of the Dirichlet-to-Neumann Map for Gratings 

In this section we will show that the scattering data from propagating modes 
determine the Dirichlct-to-Neumann map on a line x^ = T for the case of gratings. 
To do this we will first show that the traces of an appropriate family of distorted 
plane waves on x 2 — T are dense in L 2 (Q < X\ < 2ir). 

To begin we need the incoming and outgoing fundamental solutions for — A — k 2 
on R 2 in a form compatible with (1). Using Fourier series in x\ to reduce this to 
an ODE in X2, one computes that for Im{/c} > 



u + (x, k, n) = cxp(i(n + a)x 1 — ix 2 \/k 2 — (n + a) 2 ) + v + (x, k, n) 



and 



u+(x,k,n) = cxp(-ix 2 \/nn(k))(f) n (x, k) +v+(x,k,ri), 



l(-A-k 2 )-'f](x) = 




,i[(m+a)(xi-3/i)+A m (fc)|x2-!/2 



ii f(y) 



dyi)dy 2 



(3) 



where 



A m (fc) = y/k 2 - (m + a) 2 
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with the branch chosen ~ k near infinity and the cut on (— \m + a\, \m + a\). Note 
that Im{A m (fc)} > for Im{fc} > 0, and hence (-A - k 2 )- 1 maps H^(R 2 ) into 
H 2 (R 2 ). The continuous extension of A m (fc) from Im{/c} > to the real axis is 
given by 

A m (fc) = k^l-(m + a) 2 /k 2 

when (m + a) 2 < k 2 and 

\ m {k)=i^{m + a) 2 -k 2 

when (to + a) 2 > k 2 . The corresponding extension of (—A — fc 2 ) -1 to R\{±(m + 
a), to e Z} gives the outgoing fundamental solution, G + (k). For the incoming 
fundamental solution we take 

A m (fc) = -y/k 2 - (m + a) 2 , 

i.e. the branch chosen ~ — k near infinity. Substituting this in the formula for (— A— 
& 2 ) -1 / to get (— A — fc 2 ) -1 in Im{fc} < 0, and define the incoming fundamental 
solution, G-(k), by continuous extension from Im{/c} < to the real axis. Hence, 
by construction G+(fc) extends analytically to (—A — fc 2 ) -1 in Im{fc} > 0, and 
G-(k) extends analytically to (—A — fc 2 ) -1 in Im{fc} < 0. 

Now we turn to the construction of generalized distorted plane waves for L. 
Choose tp <E C°°(R) such that ip = 1 on a neighborhood of {\x2\ > T}, and the 
support of ip(x 2 ) is contained in the set where L — — A and a = 1. An (outgoing) 
generalized distorted plane wave for L is a solution of (L — k 2 )u = in D of the 
form u + = ip(x2) exp(il ■ x) + w + , where l\ = a mod 1 with < a < 1, fo = 
— \Jk 2 — (n + a) 2 , and v + is defined by limiting absorption, i.e. 

v+ = - lim (L - (k + ^e) 2 )- 1 ( J L - k 2 ){ifj(x 2 )e a - x ). 

These are generalized distorted plane waves in the sense used in §1, since the second 
component of I is not necessarily real. Ordinarily, outgoing distorted plane waves 
are defined as solutions of the form u + = exp(il ■ x) + v + where v + is outgoing. 
However, we have 

u + = ip(x 2 )e lt ' x +v + 
= e a - x + [(<P(x 2 ) - l)e ihx - lim (L — (k + ief)-\L - fc 2 )(^(x 2 )e 1 ^)]. 

Since (ip(x 2 ) — 1) exp(iZ • x) is outgoing, the term in brackets is outgoing. 
The limit defining v + will exist unless 

i) k is one of the "thresholds", k 2 = (n + a) 2 , where G±(k) are undefined, or 

ii) there is a solution to the homogeneous equation (L — k 2 )u = in D which is 
square- integrable on D n {0 < x\ < 2ir}. 

We denote the set of exceptional fc's defined by i) and ii) as S. 1 

1 Note that case ii) can occur. Choose V with compact support so that — d% 2 +V{x2), considered 
as a Schrodinger operator on IR, has a bound state, u € L 2 (R), i.e. (— <9^ 2 +V{x2))u = Eu. Then, 
taking m large enough that (to + a) 2 + E — V is strictly positive, defining e^ 1 (x) = {(m + a) 2 + 
E)((m+a) 2 +E-V)- 1 and j/> = exp(i(m+a)x 1 )u(x 2 ), we have (l/e(a;))A^+((m+a) 2 +E)ip = 0. 
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Since h = n+a for a unique n E Z, we use n and k to parametrize the generalized 
distorted plane waves, u — u(x,n,k). With these definitions we have outgoing 
distorted plane waves for all (n, fcjeZx M\S. The analytic properties of G + (k) 
discussed above carry over to the u + (x, n, fc)'s: they have analytic continuations 
to Im{&:} > which extend continuously back to R\S. This leads directly to the 
following conclusion which we state as a lemma. 

Lemma 0: For each n the set {u + (x,n, fc), k e /}, where / is an open interval in 
k 2 > (n + a) 2 determines u + (x,n,k) for k e R\S. Thus the true distorted plane 
waves determine the generalized distorted plane waves. 

The following observation is the main step in the proof. 

Lemma 1. Letting (L — (k — i0) 2 )~ 1 g denote lim e ^ _ (L — (k + ie) 2 )~ 1 g 7 the 
"incoming" solution, we have 

r 2 



/ f(xi)u + (xi,T,m, k)dx\ = 
Jo 



f e i[(m+ a )x 1 - X2 x m (k)]^ L _ k 2^ L _ (k-iO) 2 )- l (f&r)dx, (4) 

^Z5n{0<xi<27r} 

where 5r(<f>) = Jo* <t>{xi , T)dx x , f{ Xl ) e L 2 (0 < x x < 2tt). 
Proof: We have 

u+ = ^e u ' x -{L-{k + ^0) 2 )- 1 (L - k 2 )(^e hx ). 

Hence, letting D = D n {(xi,^) : < x\ < 2*7r}, for any smooth g satisfying (1) 
with bounded support in X2 

[ u+gdx= [ iPe Ux gdx- [ g(L — (k + zO) 2 ) _1 (L — k 2 )(ipe il ' x )dx. (5) 

JD JDq J D 

Since (L - (k — iO) 2 ) -1 is the adjoint of (L — (k + zO) 2 ) -1 , we have 
f g(L-(k+iO) 2 )- 1 (L-k 2 )(ijje ll - x )dx ^ f (L - {k - iQ) 2 )- l g{L-k 2 ){^e llx )dx. 

JDo JDo 

(6) 

Since L = — A on the support of ip, for any smooth h satisfying (1) 



h{L - k 2 )(ipe il - x )dx = / e ilx {2W^ • V + Aip)hdx. (7) 
Beginning with (5) and using (6) and (7), we have 

f u+gdx= [ 4>e il - x gdx- [ (L-(k- iO) 2 )" 1 ^!, - k 2 )(ipe il ' x )dx = 

JDo JD J D 

[ i)je a - x gdx+ [ (W • V + A*/0(L - (fc - iO^f^e^dx = 

J Do JDo 



/ 



e ll - x (L - k 2 )tP(L - (k - iQf^gdx 



Now approximating f(xi)Sr by g of the form above gives (4). 

With (4) we can easily prove 

Lemma 2: Assume that k € R\(5uSt) is fixed, where St is the set of k for which 
there are nontrivial solutions to Lu — k 2 u = which vanish on x 2 = T and are 
square- integrable on DC\{x 2 < T}. Then the linear span of {u + {x\,T, m,k),m e 
is dense in L 2 (0 < x\ < 2n). 

Proof: Suppose that / G L 2 (0 < ii < 2ir) is orthogonal to the span of 
{u+(xi, T, m, k), m 6 Z}. Then (4) implies 

/OO r-2-K 
( / eW'i-^Wp _ fc2)^,( L _ (fe - iO) 2 )- 1 (/<5 r )da;i)r2a:2 - 0, 
-oo J0 

for all meZ. Let 

w= de{ (L-ik-zOfy'ifST). 
Since w is incoming, we have ipw = G-(k)(L — k 2 )ipw. Moreover, we have 

(-A - k 2 )ipw = fS T - (2V?A • V + Aip)(L — (k — iOfy^Sr = 

for X2 > T. So when we represent ipw as G-(k)(L — k 2 )ipw using the analog of 
(3) for G-(k), the integrand is supported in y 2 < T. Therefore, when x 2 > T, 
\ x 2 — 1)2 1 = x 2 — J/2 on the support of the integrand in (3), and the identity above 
implies 

ipw(x) = 

for x 2 > T, i.e. w(x) = for x 2 > T. 

At this point the arguments for Case 1 and Case 2 separate. In Case 1, w is an 
incoming solution to the homogeneous problem (L — k 2 )w = on D n {x 2 < T}, 
satisfying (1) and w(x u T) = 0. Thus io=0on all of d(D n {x 2 < T}). In this 
case we have for R sufficiently large 



= / ((Lw — k 2 w)w — w(Lw — k 2 w))a(x)dx 

J £>n{0<xi <27r}n{ -_R< |x 2 1 <T} 

f27r dw dw 



ir—- Tc—dxi\ X2 =- R = Yl 2™y/k 2 - (m + a) 2 |a m | 2 + 0(e SR ). 



I 

J ° 9x2 9x2 (m+a) ^ 

(8) 

The last equality comes from the representation of ipw as G-(k)((L — k 2 )tpw), i.e. 
for x 2 < -T 



{m:(m+a) 2 </i; 2 } 

From (8) it follows that the coefficients a m of the propagating modes in w vanish, 
and w e L 2 (D n {0 < x 1 < 2tt} n {x 2 < T}). In other words w is a Dirichlet 
eigenfunction for L in D n {x 2 < T} with the periodicity condition (1). 
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In Case 2 the situation is simpler. In this case one sees immediately that (L — 
(fc — iti) 2 ) -1 (f 8t) is a eigenfunction for L on D C\ {x 2 < T}, satisfying (1), and the 
proof is complete. 

Let A(fc) denote the Dirichlct-to-Neumann operator 

. / 7 \ t du m 
A(k)h = - — on x 2 = T, 
ox 2 

where u is the outgoing solution to the boundary value problem Lu — k 2 u = in 
DC\{x 2 < T}, u = h on x 2 = T. Solutions to Lu—k 2 u = which vanish on x 2 = T 
and are square-integrable on D n {0 < x\ < 2ir} n {x 2 < T} are eigenfunctions of 
lonDfl {x 2 < T} with the periodicity condition (1). When fc 2 <E St, the set of 
eigenvalues for Lonflfl {x 2 < T}, the Dirichlct-to-Neumann map with data on 
x 2 = T is not defined. Since the line x 2 — T is chosen more or less arbitrarily, for a 
fixed k one move fc 2 out of St simply by shifting T. The set S, however, is intrinsic 
to the problem. 

If the generalized distorted plane wave u+(x, fc, m) is known for x 2 > T, then we 
know du/dx 2 on x 2 — T, and Lemma 2 has the following corollary. 

Corollary 1. The set of generalized distorted plane waves {u + (x, k,m),m € 
Z,x 2 > T} for fixed fc e R\(S T U S), determine A(fc) on x 2 = T. 

We want to recover A(fc) from the propagating modes. That follows easily at 
this point. 

Lemma 3. The scattering data from propagating modes, 

{a m (n, fc) : (m + a) 2 < fc 2 and (n + a) 2 < fc 2 } 

given for all fc € M\S, determine the distorted plane waves in x 2 > T. 
Proof: By (2) 

p -ix2\/ A m (fe) r27t 

a m {n,k) = - — e- ix ^ m+ ^v + {x u T,k,n)dx u 

it follows that a m (n, fc) is analytic in fc on the set where v+ is analytic in k. For fixed 
m and n, a m (n, fc) will be part of the scattering data from propagating modes when 
fc is sufficiently large. Thus for each m and n the scattering data from propagating 
modes determine a m (n, fc) on R\S. Thus by (2) the propagating modes determine 
the generalized distorted plane waves. 

Combining Corollary 1 and Lemma 3, we conclude the the propagating modes 
determine A(fc) for fc e \(S U S T )- 

§3. Determination of the Dirichlet-to-Neumann Map for Wave Guides 

The arguments of the preceding section apply to the wave guides with modifica- 
tions that we give here. 

Since now the unperturbed operator is — CgA, we need to replace (3) with a 
representation for the outgoing fundamental solution for — CqA. To obtain this 
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representation we separate variables and use expansion in the eigenfunctions (chosen 
to be real-valued) of the Sturm-Liouville problem 

k 2 

Cq(Xi) 

with (j) m (0,k) = 0, 4>' m (B,k) = 0. Using this basis and assuming that k is chosen 
so that H m {k) /0,meN, one checks that for / with bounded support in x 2 

u(x,k) = f] i_ / e^l..-^^^,^^,).)^^^, ( 9 ) 
^ 1 2ty/fi m (k) J[o,b]xr c o(J/i) 

is a solution to (L — k 2 )u = f when k is real. To see that this is the outgoing solution 
we will show that u(x, k) continues to a square-integrable solution when k moves 
into the upper half plane. Since the boundary conditions make d 2 jdx\ + k 2 /cq(xi) 
self-adjoint when k is real, the functions <f> m (xi,k) and /J, m (k) are analytic in k 
by Rellich's theorem. This is an elementary result here, since /x m (fe) is a simple 
eigenvalue when k is real. Thus for e > 0, if we can show that Im{/x m (A; + ie)} > 
when fj, m (k) > 0, the choice of \J fi m (k + ie) that we use here (see the definitions 
preceding Theorem 1 in §1) will make lm{^n m (k + ie) > 0. However, this follows 
immediately from the observation that d/i m (k)/dk > for k real. Thus we conclude 
that for all / for which (9) is a finite sum, u extends to a square-integrable solution 
to (L — k 2 )u = f as k moves into the upper half- plane. Thus, on the complement 
of the thresholds the operator G + (k), defined by 

G + (k)f=^2 <Pm(xi,k) / ^ m{yuk )j^-dy, 

coincides with the limit of (— CqA — k 2 !)^ 1 as Im{fc} — > 0+ on a dense set of /. 
Since an easy limiting absorption argument shows that lim e ^o + ( — CqA — fc 2 /) -1 / 
exists for / with bounded support, it follows that G + (k) is the outgoing fundamen- 
tal solution. The same construction, replacing the square roots in (9) with their 
complex conjugates, leads to the incoming fundamental solution G-(k). 

As stated in §1, distorted plane waves for the wave guide are obtained by solving 
(L — k 2 )u = with the given boundary conditions on x\ = and x\ = B for 
u = $(x, /c,m) + v+, where $ is a generalized eigenfunction for — Cq(xi)A, i.e. 
$(x,m, k) — exp(ix 2 \J Hm{k))4> m {xi, k)). Note that for this to be a true distorted 
plane wave fi m {k) should be positive. However, as in §2 we allow "generalized" 
distorted plane waves where fi m (k) < 0. As in §2 the construction of outgoing 
distorted plane waves w+(x, m, k) — <I>(x, m, k) + v' + (x, m, k) is done by limiting by 
absorption. As in §2, u + has a representation u + = ip(x 2 )& + v + where 

v+ = - lim (L-(k + ie) 2 )-\L - fc 2 )(^(x 2 )$), 

£->0 + 

with L = -c 2 A. Here the cutoff function ip e C°°(R) satisfies ip = 1 for |x 2 | > T+l 
with support contained in |x 2 | > T. As before, the functions u + do not depend on 
ip. Moreover, the exceptional set S is again the union of the thresholds and the set 
of k for which there is a nontrivial, square-integrable solution to (L — k 2 )u = in 



10 



[0, B] x R. The generalized distorted plane waves u+ have analytic continuations 
to Im{fc} > 0, and hence as in the case of gratings, {u + (x, m, k), k G /}, where I 
is an open interval in {k : fJ, m (k) > 0} determines u + (x, m, k) for k e R\S (note 
that H m {k) > for k sufficiently large for each m). In other words the generalized 
distorted plane waves are again determined by the true distorted plane waves via 
analytic continuation. 

The analog of (4) for wave guides is 

b 

f(xi)u + (xi,T 7 m,k)dxi = 



I $(x, m, k)(L - k 2 )^(L — (k — iOy^ifd^dx, 

i[0,B]xl 

and this identity shows that Corollary 1 holds for wave guides. Likewise we have 
the following analog of Lemma 3. 

Lemma 4. The scattering data from propagating modes, 

{b m (n, k) : fi n (k) >) and fi m (k) > 0} 

given for all k € R\5 I , determine the distorted plane waves in X2 > T. 
Since (2') gives, 

b m (n,k) = e - ix *^W [ ( p m { Xll k)v + {x l ,T 1 k,n)dx 1 , 
Jo 

it follows that a m (n,k) is analytic in k on the set where u+ is analytic in k, the 
proof of Lemma 3 applies here, and again conclude that the propagating modes 
determine A(fc) for k e R\(S U S T )- 



L 



§4. Reduction to the hyperbolic inverse problem 



We will begin with the wave guide problem. Consider the hyperbolic equation 

v tt = c 2 (x)Av 1 t > 

in {(xi, £2) € [0, B] x (—00, T]} with zero initial conditions, v(x, 0) =0, v t (x,0) = 0, 
and the boundary conditions 

dv 

v(0,x 2 ,t) = 0, -Q^-(B,X2,t) = 0, and v(xi,T,t) = g(xi,t). 

Let Ah denote the hyperbolic Dirichlet-to-Neumann operator corresponding to this 
initial-boundary value problem: 

dv 

A H g= 7r - for ( Xl ,t) € [0,B] x [0,oo). 
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The following theorem is a particular case of results in [B] and [KKL] (see also [El], 
[E2]). 

Theorem 2. The hyperbolic Dirichlet-to-Neumann map, Ah on x 2 = T, uniquely 
determines the sound speed c(x). 

To deduce Theorem 1 from Theorem 2 we proceed as follows. Let A(fc) be the 
elliptic Dirichlct-to-Neumann operator defined previously, for c 2 (x)A, i.e. A(k)h = 
du/dx 2 on x 2 = T, where u is the outgoing solution to the boundary value problem 

c 2 {x)Au + k 2 u = in [0, B] x (-co, T], u = h on x 2 = T 

with the zeroDirichlet boundary condition on x\ — and the zero Neumann con- 
dition on x\ = B. As we observed earlier, A(fc) is analytic in k off the discrete set 
5U5t- Hence, using the Fourier-Laplace transform in t ,we can recover Ah from 
A(fc), given for k — e < k < k + e. Since we showed in §3 that the propagating 
modes determine A(fc), this completes the proof of Theorem 1 for wave guides. 

Since we have also shown for gratings that A(k) for k e R\(SUSt) is determined 
by scattering data from propagating modes, the only change in the argument needed 
to prove Theorem 1 for gratings is in the citation of results on the hyperbolic 
Dirichlet-to-Neumann map. Here the hyperbolic Dirichlet-to-Ncumann operator, 
Ah, is defined by 

= q — on x 2 = T, 
ax 2 

where v is the solution to 

v tt = Lv in D n {x 2 < T} x {0 < t < oo} 

satisfying the periodicity condition (1), the initial-boundary conditions v(x,0) — 
v t (x, 0) = and v(x\,T, t) = g. In this setting the uniqueness results of [B], [KKL] 
and [El, 2] imply that that Ah given on x 2 — T determine both the coefficients 
of L, i.e. the pcrmitivity e(x), and the domain D. This completes the proof of 
Theorem 1. 
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